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ABSTRACT 

We present here the results of cross-correlating the X-ray background 
measured by Ginga in the 2-10 keV band with several catalogues of 
extragalactic objects. Positive signals with an amplitude of a few per cent 
have been found for some catalogues implying that some fraction of the 
X-ray background is produced either by the class of catalogued sources 
or by other classes spatially related to them. Detailed X-ray background 
simulations have been used to assess the significance of the results and, 
for the first time, the full angular shape of the cross-correlation. The 
inferred X-ray volume emissivity in the local Universe, jo: has been 
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estimated for two galaxy catafogues (UGC and IRAS) for which the 
cross-correlation is highly significant. We obtain jo = (0.74 ± 0.07) for 
UGC and jo = (1.15 ±0.10) for IRAS, in units of lO^H erg s'^ Mpc-^. 
Extrapolating this result back to 2; ~ 1 — 4 leads to the conclusion that 

< 10 — 30 per cent of the X-ray background could be produced by a non- 
evolving population of galaxies. These values are shown to be consistent 
with upper limits on the Auto Correlation Function derived here. 



2 



1 INTRODUCTION 



The puzzle of the origin of the X-ray background (XRB) still remains unsolved, al- 
though it is becoming apparent that no single class of 'miraculous' sources is able to 
satisfy simultaneously all the observational constrains accumulated over the last 30 
years (see Fabian and Barcons 1992 for a review). 

The fraction of the XRB directly resolved into sources in the 'soft' band (< 
2 keV) is about 50 per cent in the deepest ROSAT surveys (Hasinger et al. 1993, 
Brand uardi-Raymont et al. 1994). The brightest of these sources (i.e., those with 0.5- 
2 keV fluxes > 10~^^ erg cm~^ s~^) are expected to be quasars with z < 2.5 (based on 
the optical identification work of Shanks et al. 1991). Optical identification programs 
are being carried out in order to establish the nature of the sources at and below that 
flux. On the other hand, fluctuation studies of these Deep Rosat fields (Hasinger et 
al. 1993, Barcons et al. 1994) reveal the imprint of ~ 1000 sources deg~^, 'resolving' 
a fraction of ~ 70 per cent of the soft XRB. The nature of the sources at fluxes 
~ 10"-^^ ergcm"^ s~-^ will not be known until X-ray telescopes more sensitive than 
Rosat become available. 

The lack of imaging devices in the 'hard' band (2-10 keV) has made the directly 
resolved fraction of the XRB in that band much smaller, only a few per cent (Piccinotti 
et al. 1982). Integrating and extrapolating the obtained luminosity functions, clusters 
of galaxies would contribute ~ 4 per cent of it and AGN (mostly Seyfert galaxies) ~ 20 
per cent. Although soon ASCA, and SPECTRUM-X-T in the future, are expected 
to give some insight into this problem, useful information is available contained in the 
angular distribution of the intensity of the XRB (fluctuations and anisotropies). 

Fluctuation studies using Ginga data (Butcher et al. 1994) have shown that 
the log — log S curve is consistent with an euclidean shape (sources per unit flux 
and unit solid angle n{S) oc S~'^'^) down to ~ 10~^^ ergcm"^ s~^ (2-10 keV). Recent 
XRB anisotropies and excess fluctuations studies have shown that only 60-70 per cent 
of the 'hard' band XRB can be produced by sources clustering on scales ~ 6/i~^ Mpc 
{Hq = lOO/ikms"-^ Mpc~^) (see Danese et al. 1992 for a recent review). This implies 
that the obvious solution of just 'adding' more QSOs until the XRB saturates does 
not trivially work. 

A complementary way of facing the problem of the origin of the XRB is to 
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search for cross-correlations of the XRB intensities with positions of known sources in 
existing catalogues. 

We define the cross-correlation function (CCF) as 



_ m ■ 5I)e 
^^^^^^ - {N){I) ' 

5N being the fluctuation in projected galaxy density and 51 the XRB deflection, (A^) 
and (/) being the mean projected galaxy density and XRB intensity, respectively (see 
Sections 2 and 3 for definitions). {)q means an average among all those pairs distant 
an angle d. 

Jahoda et al. (1991) obtained a value of 3 x 10~^ for the zero-lag cross- 
correlation function (CCF) of the 2-10 keV XRB intensities (as observed by HEAO-1) 
with ESO and UGC galaxies. Jahoda et al. (1992), using this value, concluded that a 
sustantial fraction of the 'hard' XRB could be produced by a non-evolving population 
of X-ray sources. 

Lahav et al. (1993) performed a more sophisticated but preliminary analysis 
of Ginga and HEAO-1 data at zero-lag. They arrived at the conclusion that ~30-50 
per cent of the 'hard' XRB could be produced by local non-evolving X-ray sources. 
Their analysis (unlike that by Jahoda et al. 1992) took into account the clustering of 
the sources, which, as shown there and in Section 4, can significantly reduce the values 
of the emissivity deduced from the observed CCFs. 

Miyaji et al. (1994) obtained jo = (0.9 ± 0.2) x 10^^ h erg s-^ Mpc-^ using 
the zero-lag CCF of the whole sky HEAO-1 XRB observations with two catalogues of 
IRAS galaxies. They also used a sample of AGN detected by HEAO 1 (Grossan 1992) 
to study the spatial correlation function and the selection function of sources detected 
both in X-rays and infrared wavelengths. 

In the present paper we have analyzed the CCF obtained from a set of Ginga 
scan data (see Section 2.1) with nine catalogues of extragalactic sources (listed in 
Section 2.2). We have used, for the first time, the full angular shape of the CCF up 
to 5 degrees to gather information about possible sources of the 2-10 keV XRB and 
their relation to the sources in the catalogues. We have also taken into account the 
limits imposed by the observed auto-correlation function (ACF) of the XRB over the 
same region of the sky. 



In Section 3 we show how the CCF of the XRB with the galaxy catalogues is 
obtained. We also discuss the significance of the signals seen, basing our arguments on 
detailed simulations of the X-ray sky and the galaxy catalogues in the region under 
study (Section 3.2). 

Section 4 is devoted to the theoretical framework in which the CCF is modelled. 
Wc give the expression of the expected CCF as a function of the properties of the 
sources (selection functions, spatial correlation functions and emissivities) and the 
collimator profile ('beam'). 

Our results are used in Section 5 to estimate the volume emissivity of X-ray 
sources related to a local population of galaxies. The implications of these results are 
then addressed in Section 6. 

Finally, in Section 7 we summarize our results. 

2 THE DATA 

2.1 The X-ray data 

Our X-ray sample consists of 5 strips over the North Galactic Pole observed in scan 
mode with the Large Area Counter (LAC) on board of Ginga. Every strip was scanned 
several times by the satellite, accumulating counts in 16 s time intervals, which cor- 
respond to ~ 0.2 degrees angular separation. As consecutive scans in the same strip 
do not exactly overlap, we actually have measurements of the brightness of the X-ray 
sky with relative separations smaller than that distance. However, as the collimator 
profile is only well determined in 0.1 deg steps, we have binned the observed counts in 
0.2 deg wide bins. This binning also helps in improving the signal-to-noise ratio. 

A more thorough description of the data and their reduction process is given 
in Carrera et al. (1993), our data being a subset of those used in that work. We 
only stress here that the energy band used {Ginga LAC Pulse Heigth Analyzer -PHA- 
channels 8-20, corresponding roughly to 4-12 keV), and the high galactic latitude 
of the sample (| & |> 30 deg), lead to a small contamination by diffuse gas from our 
galaxy. Furthermore, the extension of this region is small and the angle between it and 
the dipole direction is 86 degrees, making insignificant any dipole correction because 
of the Compton-Getting effect {e.g., Boldt 1987). Given that our sample is at a 
supergalactic latitude ~ 40 degrees, the supergalactic plane do not pose any problems 
of fair sampling on it. 
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As the region of the sky scanned is also away from known hard X-ray sources 
(those of Piccinotti et al. 1982), we have a sample of extragalactic diffuse XRB 
deflections or fluctuations 51 = I — (/) , (/) being the mean XRB intensity. 

2.2 The catalogue samples 

A total of 7 galaxy samples have been used in this work, together with the Abell 
catalogue of clusters of galaxies. The sources from each catalogue with Right Ascension 
(a) e [210, 270] degrees and Declination (5) G [30, 70] degrees have been selected and 
their positions used to define a set of fluctuations in projected galaxy density (SN) 
in the following way: at the same positions and with the same collimator orientation 
as in the real XRB observations, a series of galaxy number density (AT) 'observations' 
have been performed. Then they are binned in 0.2 deg wide bins (again in the same 
way as in the X-ray observations) and their mean value {{N)) subtracted to give the 
fluctuations 5N. 

The number of sources used from each catalogue and the value of (A^) for each 
one of them (in sources beam~^) are given in Table 1. 
Briefly the catalogues are as follows: 

- IRAS: the Meurs and Harmon (1989) catalogue that comprises infrared galaxies 
selected (using purely colour criteria) from the IRAS Point Source Catalogue 
(PSC) with 60 /im fluxes > 0.7 Jy. This catalogue was also used by Miyaji et 
al. (1994). 

- UGC: the Uppsala General Catalogue (Nilson 1973) of optical galaxies, selecting 
those whose major angular diameters are greater than 1 arcmin. 

- IRASQDOT: IRAS QDOT catalogue taking one in six sources out of a galaxy 
subsample of the PSC (Rowan-Robinson et al. 1990). 97 per cent of the 2163 
sources in this catalogue (60 //m fluxes > 0.6 Jy and | & |> 10°) have been 
spectroscopically identifled. 

- CfA: Centre for Astrophysics redshift catalogue of galaxies (Huchra et al. 1983), 
which happens to partly overlap with the region of the sky observed in X-rays. 

- Abell: positions of Abell (1958) clusters. The richnesses of the selected clusters 
are between and 4, most of them having Richness class or 1. 

- IRASZ: Strauss et al. (1990) sample of IRAS PSC galaxies complete down to 
60 //m fliixes > 1.936 Jy, all with | b \> 5°. Redshifts are available for these 
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sources. There are two known AGNs in this catalogue within our sample region; 
including or excluding them does not make any appreciable difference to the 
results discussed below (neither does it to those of IRASZS or IRASZSVOL 
below) . 

- IRASZS: a selection among galaxies in the IRASZ sample, choosing those with 
500 < c^LocaiGroup < 8000 km s~^. This is done (Yahil et al. 1991) to avoid 
both local galaxies with possible peculiar velocities and very distant ones, for 
which the selection function (probability of a source being detected; see Section 
4) is very small. This catalogue was also used by Miyaji et al. (1994). 

- IRASZSVOL: this comprises the same sources as in IRASZS, but we weighted 
them by the inverse of the selection function at their distances. The idea behind 
this weighting is to 'compensate' for the sources missing in the sample (the 
farther away, the more likely to be missing, the smaller the selection function 
and the larger the weight). 

3 THE OBSERVED CROSS-CORRELATION FUNCTIONS 
3.1 Measurement of the Cross Correlation functions 

In Section 1 we defined the CCF as 

_ m • 5I)e 
^^^^^^ - {N){I) 

To evaluate (/) we have used the spectrum of the XRB obtained by Marshall 
et al. (1980). It leads to (7) = 5.7 x 10"^ ergcm-^ s"^ sr-^ in the 2-10 keV band, 
which corresponds to ~ 8ct s~^ beam~^ in the LAC PHA channels 8 to 20. 

The values of VFxg(O) in 1 deg bins derived from the different samples are noted 
in Table 1. The UGC-XRB CCF value is larger than the one obtained by Jahoda et al. 
(1991) (W^xg(O) - (2 - 5) X 10-3). Our CCF value for the IRAS and IRASZS are also 
larger than their equivalents in Miyaji et al. (1994). This is hardly surprising, since 
the CCF is heavily dependent on instrument characteristics (/ is measured in counts 
beam"-^ and N in galaxies beam"-^). For example, the Ginga LAC probes deeper in 
source counts than HEAO-1 A2 due to its smaller beam size (2° x 1° and 3° x 1.5° 
respectively) . The different beam sizes also have a geometrical effect that contributes 
to make the signal different (see Section 4). 



We have plotted Wxg(^) in 1 deg wide bins for the UGC and IRAS samples 
in Figs, la and 16, respectively (filled dots). The error bars are from the standard 
deviation in each bin. As we will show in Section 4, CCF signals (such as those seen 
in Figs, la and lb) are expected due to the finite size of the Ginga LAC collimator 
(~ 1° X 2°) and/or to the spatial clustering of the X-ray sources with the galaxies in 
the samples (see also Lahav et al. 1993). 

3.2 Significance of the Cross Correlation Signal 

We have simulated both catalogue samples (simply by reshuffling at random the po- 
sitions of the objects in them) and X-ray fluctuations to assess the signiflcance of 
our results. By obtaining the CCFs of these random samples, we can test if purely 
geometrical effects or random positioning could give rise to the observed CCFs. 

The simulated X-ray samples have been obtained by uniformly distributing X- 
ray sources following the log N — log S curve obtained by Butcher et al. (1994) (see 
also Carrera et al. 1993 for a complete description of the simulation process). 

All these simulated populations (both X-ray and catalogue) have then been 
folded through the Ginga LAC collimator and the corresponding deflections have been 
obtained as described in Section 2.2. 

The values in the column labelled VFxg,sim(0) of Table 1 are obtained using 100 
simulations of every galaxy sample and 100 simulated XRB samples: the quantities 
quoted are 2cr upper limits (they are the higher of the two values that encompass 95 
per cent of the simulated VFxg(O) for each sample). The CCFs from the samples UGC, 
IRAS and IRASZS are well above those upper limits, showing that they are >> 2cr 
significant, and that they do not come from some chance positioning of the sources in 
the samples or the X-ray data, nor from some geometrical effect due to the collimator 
or the way the CCFs have been obtained. The signals for the CfA and IRASZ samples 
are just about at the 2cr significance level. 

These differences between the CCFs from different samples reflect their dif- 
ferent densities, flux limits and clustering properties, and, of course, their different 
percentages of/clustering with X-ray sources giving rise to the XRB. 

In Figs, la (UGC) and Ih (IRAS) the medians of the simulated CCFs 
(W/xg,sim(^)) are shown as solid lines, as well as values encompassing 68 per cent (dashed 
lines) and 95 per cent (dotted lines) of the simulations. 
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To test whether most of the signal is due to a few sources just unresolved in 
X-rays, we also calculated the CCF (both real and simulated) excluding the five per 
cent higher X-ray deflections. Both UGC and IRAS CCFs were still above 97 per 
cent of the simulations, but IRASZS (and CfA and IRASZ) was well below that limit. 
This, as said earlier, indicates that the CCF signal found for IRASZS was due to a few 
'bright' X-ray sources just unresolved, while the UGC and IRAS signals come from 
the sources in those samples as a population. Therefore, we pay special attention to 
the understanding of these signals. 

The simulated CCFs have also been used to estimate the error bars for the 
observed CCFs at each angular separation (from the 68 per cent limits for each simu- 
lation). These are shown in Figs. 2a and 2b (see Section 5). 

Having checked and established the significance of the observed CCFs, their 
relationship to the properties of the (catalogue and X-ray) sources will be outlined in 
next section. 
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4 THE THEORETICAL CROSS-CORRELATION FUNCTION 

In order to model the CCF on the basis of the relationship between the X-ray sources 
responsible for the observed fluctuations in the X-ray sky and the catalogued opti- 
cal/IR sources, let us consider two (distinct) source populations. The first one would 
consist of galaxies with selection function P{r) (probability of a source being de- 
tected at a distance r from the observer) and mean spatial density {fig). Our second 
population would consist of X-ray sources with local volume emissivity jo (emitted 
X-ray power per unit volume) and clustered with the above galaxies, with a spatial 
cross-correlation function Cxg(rg — i"x) (probability in excess of a purely poissonian 
distribution that a galaxy and an X-ray source lie at a distance | rg — Fx |). We will 
assume, as usual, a form ^xg(?') = {f/fo)~^ ^or Cxg{f)- Cxg{r) — (or equivalently. 
To — 0) would mean that these two populations are independently distributed in the 
sky. 

Some of the X-ray sources may be actually present in the catalogue. In this 
case X-ray and catalogue observations in different directions (given by the unit vectors 
Hx and Hg respectively) would be correlated when viewd through a collimator with 
profile G(nx — Hg). If we define 



then the correlation term without taking into account the possible clustering of the 
catalogued sources (poissonian term) would be 



ri{n^ - ng) = (/) • (N) ■ W^^in^ - ng) 



(2) 




dr r 



,2 J^g 




(3) 



where jxg jo) is the local X-ray volume emissivity of the galaxies that are also X- 
ray emitters (we have assumed that jxg is constant within the galaxy sample depth) 
and Pyi{r) their selection function. Defining the effective catalogue depth as 




(4) 



and the collimator autocorrelation function (ACF) as 
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then 



?7pois(nx - ng) = ^ixgi?*r7coi(nx - ng) (6) 
If we now allow for clustering among the two populations, the CCF is 



?7cius(nx - ng) = J dV^' j dVg' P{r^'){n^') G{n^> - nx)G(ng/ - ng)^(rx' - rg/) 

(7) 

which, under the small angle approximation and using Eq. 5, leads to 

^7cius(nx-ng) = ^jo{ng)H^r^ J drr^~^P{r) J d'^Vtn r7coi(n-nx) | n-ng |^"^ (8) 

where again we have assumed a constant jo over the volume of the galaxy sample and 
= r(l/2)r((7 - l)/2)/r(7/2), r being the gamma function. 

The total CCF would then be the sum of these two terms (Eqs. 6 and 8). If 
clustering is present, but its contribution neglected, the volume emissivity is grossly 
overestimated (by as much as a factor of 6 in the relevant case). Eq. 8 generalizes 
Eqs. 2 and 3 in Lahav et al. (1993) paper for the non-zero lag non-square collimator 
case. Note that the emissivity {jq in this paper) was denoted as px in that paper. 

For the UGC and IRAS samples the selections functions are well determined 
(Hudson and Lynden-Bell 1991 and Yahil et al. 1991 respectively). 

Miyaji et al. (1994) constructed the spatial correlation function of a sample 
of HEAO 1 X-ray sources with the IRAS catalogue, obtaining 7 = 1.8 and Tq = 
4h~^ Mpc. They also derived the selection function for the X-ray sources that were 
also in the IRAS 0.7 Jy sample (of which our IRAS sample is a subset) Px{f), obtaining 
a value of i?* = 65 Mpc. These are the values of 7, tq and i?* that we have used 
for our IRAS sample. 

For the UGC sample we have assumed Px{r) = P{r) (the selection function of 
the X-ray emitting galaxies in the UGC catalogue is the same as the global one) for 
simplicity, since its real value is not known. Furthermore, we have used the spatial 
correlation function of the optical galaxies as the cross-correlation function of them 
with the X-ray sources, i.e., 7 = 1.8 and ro — 5h^^ Mpc. 
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Finally, in order to minimize the number of free parameters, we have assumed 
jxg = jo (the X-ray sources are all galaxies, and they are clustered like the galaxies of 
the catalogue). The effects of this assumption are discussed in Section 6. 

Taking into account all the above assumptions, from Eqs. 6 and 8 we obtain 

?7(nx - ng) = ^ ( i?*?7coi(nx - ng) 



+ {n^)H^rl / d^r r^-'^P^r) / d'iln ?7coi(n - nx) | n - ng l^-^ j 
The value of (ng) is easily calculated taking into account 

{N) = J d^r P{r)G{n){ng) = neff J dr r^P{r) ■ (ng) (10) 

where ^e// = / d'^Q^ G'(n) is the effective solid angle of the collimator (in our case 
5.7 X 10~^ sr), and the values of (A^) are listed in Table 1. 

The collimator profile of the Ginga LAC is well known down to scales ~ 0.1 
degrees, but unfortunately the geometry of the X-ray observations is quite complex, 
involving many different relative positions and orientations of the collimator. We have 
then obtained r]co\{Q) from the set of simulations of X-ray sources previously used to 
estimate the significance of the observed CCFs, and using the relation {51 ■ 51) q = 
(S*^) ■ rico\{d) where (S*^) is the mean of the square of the X-ray flux S as calculated 
from the log — log S curve used in the simulations (see Section 3.2). In this way we 
obtain a collimator ACF 'tailored' to our needs, including all the geometrical effects 
of our X-ray data. 

Finally, the 'zero' level of the CCF (i.e., the cross-correlation expected in the 
absence of any true cosmic signal) is not zero because the finite size of the sample 
introduces a negative baseline in the correlation functions (Kondo 1990). This zero 
level (= — Wxg(O) / <Number of pairs at zero-lag>) has been calculated and subtracted 
from the observed CCF before fitting. 
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5 RESULTS: THE VOLUME EMISSIVITY 



Although our angular bin size (0.2 deg) provides us with many bins, they are hardly 
independent (due to the collimator extension and to the fact that we are dealing with 
the product of pairs of observations each of which appear in many pairs). Therefore, 
we have only kept one free parameter (Jq). Its value has been obtained by fitting 
(using minimum least squares) the model given in Eq. 9 to the observed CCF binned 
in 0.2 degree bins. This has been done for two difi'erent angular ranges (0.2 — > 2 deg 
and 0.2 5 deg) and for several values of -Rmin, which corresponds to the distance 
beyond which the selection function drops below one (P(r < -Rmin) = 1), ^-c, galaxies 
begin to be lost or undetected. The results are shown in Table 2. 

There is a significant difference between fits to the different angular ranges: 
this is because r7cius(^) is significantly different from zero between 2 and 5 degrees 
and then the fitting over the larger range is bound to require a lower normalization 
(— > lower jo). This is shown in Figs. 2a (UGC) and 2b (IRAS) where the models 
marked with an asterisk in Table 2 are plotted along with the observed CCFs: the 
solid lines correspond to the total CCFs, while dashed lines are the clustering terms 
and dotted lines the poissonian ones. The error bars on the observed CCFs are from 
the simulations of Section 3.2. 

The peaks visible in Figs. 2a and 2b at about 1.6 degrees in the observed CCFs 
are very probably due to an X-ray source just unresolved that happens to fall between 
two of our strips. If the simulations explained above are repeated with 0.2 deg bins, 
that peak is also seen in 5 out of a 100 simulations. A similar effect (due to the 
overlapping of the long sides of the collimator in every two adjacent strips) is seen in 
the poisson CCFs. 

From the results in Table 2, the UGC sample yields an X-ray volume emissivity 
jo = (0.74 ± 0.08 ± 0.13) X 10^^ h erg s"^ Mpc-^ while the IRAS one gives jo = 
(1.15 ± 0.10 ± 0.19) X 10^^ h erg s~-^ Mpc~^. These values are in good agreement with 
the zero-lag results obtained by Lahav et al. (1993). The errors given are the 1 and 
2(7 confidence levels obtained with bootstrap simulations: for each angular bin in Wxg 
we have extracted at random (and replaced) as many pairs (SI, SN) as in the real data 
and calculated {5N ■ 51)0/ {N) {I). This was repeated 1000 times and the resulting 
CCF was fitted to Eq. 9 to obtain jo; the quoted errors enclose 68.3 per cent (1 sigma) 
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and 95.4 per cent (2 sigma) of the simulations. They also span the range of values of 
jo in Table 2 obtained for each sample, reflecting thus not only the statistical errors 
in our fits, but also uncertainties in our model parameters (i^min? 7? '^o)- 

If we had ignored the clustering term (as did Jahoda et al. 1991), we would 
have obtained Jo ~ 8 x 10^^ h erg Mpc~^ for UGC, showing the importance of 
taking it into account. 

The value found by Miyaji et al. (1994) (jo = (0.9 ± 0.2) x 

10^^ h erg s~^ Mpc~^) is in good agreement with our IRAS result (they overlap within 
1 sigma), taking into account the different IRAS sample (our 0.7 Jy versus their 2 Jy 
with a local cutoff), the different instrument {Ginga LAC versus HEAO-1 A2), the 
slightly different observation bands (4-12 keV versus 2-10 keV) and the different tech- 
niques (full angular and collimator resolution versus zero-lag square collimator) used. 

6 IMPLICATIONS FOR THE X-RAY BACKGROUND 

In this section we study the fraction of the XRB contributed by the sources whose 
emissivity we have just obtained from their CCF with optical and IR galaxies. 

The intensity received from objects with local volume emissivity jo distributed 
up to redshift -Zmax is 

/(< Zma^) = %^^jo/(^max,fio,P,Q!) (11) 

where / is the effective look-back factor (Boldt 1982, Lahav 1992) which depends on 
the cosmology (r2o=l), on the X-ray spectrum of the sources (we assume a power law 
with energy index a = 0.7 which is consistent with the XRB fluctuations studied by 
Butcher et al. 1994), on the evolution of the emissivity (comoving j{z) = jo{l + z^) 
and, of course, on the maximum redshift of integration 2;max- Under this assumptions, 
we get 

/(< W) = ^^Jo dz (1 + ^oz)-'/\l + ^)-2-«+^' (12) 

47r Ho Jo 

The fraction of the XRB contributed by those objects is then F = /(< 
Zma.x)/{I)- In Figs. 3a and 3b we show contour plots for different values of F (for 
jo = 0.74 and 1.15 x 10^^ h erg s"^ Mpc-^, UGC and IRAS best flts respectively) in 
the (-2max,p) space. We can see that between 10-20 per cent of the XRB could be 
produced if we extrapolate the obtained emissivity to ^max = Ij and 20-30 per cent 
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if -^max = 4, in the absence of evolution (p = 0). Some moderate positive evolution 
would lead to much higher values (but see below). Note that Eq. 12 only depends on 
—a +p, not on their individual values, so if a different spectral slope a' is considered, 
the corresponding F would be given at (-2max,p' ~ p — a' + 0.7). 

The population of sources whose emissivity has been obtained would also pro- 
duce some auto-correlation signal on their X-ray intensities (ACF). Any estimation of 
their contribution to the XRB must be consistent with the current upper limits on the 
ACF. We have derived (Fig. 4) the non poissonian ACF of the X-ray sample in the 
same way as Carrera et al. (1993) (i.e., by subtracting from the observed ACF that 
expected from a purely poissonian distribution of sources, obtained by simulations). 
Again, only upper limits to the cosmic signal are obtained, a null value being the best 
estimate of the non poissonian ACF. 

If we assume an evolution for the spatial correlation function of ^(r, 2) = (1 + 
^-j-3-e^^^^^^-7^ then under the small angle approximation, the non poissonian ACF 
is given by 



W^^a;(ni - na) = (i^^/) ) ^^^^ I ^^^"^ - ni) | n - na |^"^ 

(13) 

dz (1 + Qoz)-^/^{l + zy^-'-^'^+^Pdl^^iz). 



f 

Jo 



'0 

Or, taking into account Eq. 11 



Wa:x{ni - 112) = —7^2 — 7^ \foHy / d'^Qn ?7coi(n - ni) I n - n2 \^~^ 



f 

Jo 



C ^lffP{Zma^,no,p,a) 

dz (1 + noz)-^/^{i + 2)-^-^-2"+2^dJi"^(^), 

(14) 

where dA{z) is the angular distance. In this expression, for a given pair (^max,p), the 
only free parameter is F, the fraction of the XRB produced by the objects whose ACF 
we are studying. 

To get a quantitative idea of the constraints that this ACF places on the possible 
values of p and 2;max5 we have x^~fitted the observed non-poissonian ACF (Fig. 4) 
to Eq. 14 over the range 0.2 — 5 deg, with F as a the only free parameter. The best 
fit (under the current assumptions on the clustering properties of the sources) always 
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corresponds to F = 0; for each pair (2;max)P) the maximum allowed value of F (from 
Ax^ = 4 or 2cr) is plotted in Figs. 5a and 56 for comoving and stable clustering 
(e = —1.2 and e = 0, respectively) and for UGC (7 = 1.8 and tq = 5h~^ Mpc). Figs. 
5c and 5c? correspond (respectively) to the same evolutions of the clustering for IRAS 
(7 = 1.8 and ro = 4/i~^ Mpc). Combining Figs. 3 and 5 we can see that the ACF 
upper limits do not allow contributions to the XRB greater than ~ 40 — 50 per cent, 
even with evolution. 

In Section 4 we assumed that jxg = jo- If we keep both quantities separate, jo 
(the total emissivity) is replaced by jxg (the emissivity of the sources in the catalogue) 
as the free quantity to fit, and ro by ro' = ro{jo/ jxg)^^'^ > tq. Hence, the net eff'ect of 
only a fraction of the galaxies being X-ray emitters would be to enhance the importance 
of the clustering term with respect to the poisson term (as expected). As the clustering 
term is already dominant for the UGC sample, reducing jxg/jo even down to 25 per 
cent does not have any effect on the derived total Jq. However, the poisson term is 
dominant for the IRAS sample and as we decrease Jxg/jo the clustering term takes over. 
This, in turn, increases the total emissivity (by ~50 per cent when jxg/jo = 0.25), but 
the ACF constraints keep F <20 per cent. 

The difference between the X-ray emissivities obtained for the UGC and IRAS 
samples is likely to reflect an intrinsic difference in properties of their constituent 
galaxies. For example, the UGC catalogue selects galaxies by their diameter, without 
distinction of type, while the IRAS catalogue is a far-infrared selected sample, with a 
high proportion of star forming and dust rich galaxies, hence mostly spirals. 

We know that sources producing the Ginga fluctuations have a spectrum with 
energy index ~ 0.7 (Butcher et al. 1994), and these are precisely the sources whose 
emissivity has been obtained here. With only a moderate contribution to the XRB, 
the residual spectrum of the XRB would be much flatter than the presently observed 
one (energy index ~ 0.4). This will have wide implications for the sources of this 
residual XRB, which will need to be highly absorbed and/or reflection dominated. 
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7 CONCLUSIONS 

X-ray observations of the 2-10 keV XRB in scan mode performed with the Ginga 
sateUite have been cross-correlated with several catalogues of different objects (mainly 
galaxies). The significance level of the resulting CCF signals have been established 
using simulations of both X-ray and catalogue samples. 

We have concluded that there is a strong cross-correlation between the hard 
XRB and nearby galaxies (especially those in UGC and the IRAS catalogue or Meurs 
& Harmon 1989), finding W^g{0) ~ 1.2 - 1.4 x 10"^. 

Expressions have been derived for the expected CCF at non-zero lag when 
the population of X-ray sources dominating the Ginga fluctuations are associated or 
strongly clustered with the catalogued galaxies. In the specific case of the X-ray 
sources being present in the galaxy catalogue, least-squares fitting has been used 
to obtain the X-ray volume emissivity jo necessary to produce the observed signal, 
obtaining jo ~ 0.7 - 1.2 x 10^^ h erg s"^ Mpc"^. 

Extrapolating this emissivity back to 2; = 1 in the absence of evolution, 10-20 
per cent of the XRB would be produced (<30 per cent if the contribution from the 
non-evolving sources is integrated up to 2; = 4). We have shown that these fractions 
are consistent with the absence of any ACF signal in the XRB. 

Higher fractions, close to saturation, would be achieved even with very moderate 
evolution. However, these high fractions are not permitted by the upper limits on the 
ACF. Furthermore, the spectrum of the summed contribution of these sources would 
be very different from the observed XRB spectrum, unless rapid spectral evolution 
also takes place. 

Direct observations with higher angular resolution instruments (ASCA, 
SPECTRUM-X-F, XMM) would be very useful to decide on whether the emissivi- 
ties obtained here are actually from the sources in the catalogues, or come from X-ray 
sources clustered with them. They would also provide spectra for fainter sources, 
helping to decipher the riddle of the hard X-ray background. 
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FIGURE CAPTIONS 

Figure 1. Wxg(6') in 1 degree bins (filled dots): the error bars are extracted from the 
standard deviation in each bin. Also shown are the median of the simulations (solid 
line) and the limits encompassing 68 (dashed lines) and 95 (dotted lines) per cent of 
the simulations, a is for the UGC sample and h for the IRAS sample (see text). 
Figure 2. Wxg(^) in 0.2 degree bins (filled dots): the error bars are extracted from the 
simulations. Also shown are the best total fit (solid line) and the clustering (dashed 
line) and poissonian terms (dotted line) . a is for the UGC sample and h for the IRAS 
sample. The peaks at about 1.6 degrees present in both observed CCFs are due to a 
just unresolved X-ray source between observed sky strips (see text). 
Figure 3. Contour plots (in the (-Zmax,^) space) of the fraction of the X-ray back- 
ground F produced by a population with local emissivity jq: levels of F=0.1 (solid 
line), 0.3 (dashed line), 0.5 (dot-dashed line), 0.7 (dotted line) and 0.9 (dash-dot- 
dot-dot line) are shown, a is for jq = 0.74 (UGC) and h for jo = 1.15 (IRAS), both in 
units of 10"^^ h erg s~^ Mpc~^. 

Figure 4. Non-poissonian ACF of the XRB: the error bars {la) are extracted from 
simulations (see text). 

Figure 5. Contour plots of the maximum F allowed by the ACF. Levels of F=0.1 
(solid line), 0.3 (dashed line), 0.5 (dot-dashed line), 0.7 (dotted line) and 0.9 (dash- 
dot-dot-dot line) are shown, a is for comoving clustering and UGC, h for stable 
clustering and UGC, c is for comoving clustering and IRAS and d for stable clustering 
and IRAS. 
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Table 1. Cross-correlation functions or several samples with the XRB (see text). 



Name 


N 


(N) 


xg \^ ) 


W^xs simfO) 
xg,sim v"-"/ 


UGC 


951 


1.24 


0.0116±0.0015 


<0.008 


IRAS 


351 


0.48 


0.0142±0.0016 


<0.010 


IRASQO 


534 


0.74 


0.0063±0.0013 


<0.009 


CfA 


196 


0.32 


0.0190±0.0030 


<0.018 


Abell 


210 


0.26 


0.0100±0.0020 


<0.017 


IRASZ 


82 


0.10 


0.0240±0.0050 


<0.026 


IRASZS 


45 


0.06 


0.0460±0.0080 


<0.026 


IRASZSVOL 


45 


1.97 


0.0060±0.0060 


<0.069 



Table 2. Results of the fitting for different parameters (see text). 



Sample 




7 


Range fitted 


-f^min 


io (±la ± 2(7) 




{h-^ Mpc) 




(degrees) 


(/i-i Mpc) (10^9 /i erg s-^Mpc-^) 


UGC 


5.00 


1.80 


0.2,2.0 


11.1 


0.82 


UGC 


5.00 


1.80 


0.2,5.0 


11.1 


0.65 


UGC 


5.00 


1.80 


0.2,2.0 


13.3 


0.92 


UGC* 


5.00 


1.80 


0.2,5.0 


13.3 


0.74 (±0.07 ±0.13) 


IRAS 


4.00 


1.80 


0.2,2.0 


1.0 


1.22 


IRAS 


4.00 


1.80 


0.2,2.0 


5.0 


1.23 


IRAS 


4.00 


1.80 


0.2,2.0 


10.0 


1.24 


IRAS 


4.00 


1.80 


0.2,5.0 


1.0 


1.14 


IRAS* 


4.00 


1.80 


0.2,5.0 


5.0 


1.15 (±0.10 ±0.20) 


IRAS 


4.00 


1.80 


0.2,5.0 


10.0 


1.17 


IRAS 


4.00 


1.65 


0.2,2.0 


5.0 


1.34 


IRAS 


4.00 


1.65 


0.2,5.0 


5.0 


1.23 


IRAS 


3.75 


1.65 


0.2,2.0 


5.0 


1.40 


IRAS 


3.75 


1.65 


0.2,5.0 


5.0 


1.30 
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